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We study a closed model of the universe filled with viscous fluid and quintessence 

O 

matter components in a Brans-Dicke type cosmological model. The dynamical equa- 
tions imply that the universe may look like an accelerated flat Friedmann-Robertson- 



X 



Walker universe at low redshift. We consider here dissipative processes which follow 
a causal thermodynamics. The theory is applied to viscous fluid inflation, where 
accepted values for the total entropy in the observable universe is obtained. 



I. INTRODUCTION 

Recent observation of the Hubble diagram for supernovae la indicates that the expansion 
of the universe is acceleranting at the present epoch . This acceleration is attributed 
to a dark energy residing in space itself, which also balances the kinetic energy of expansion 
so as to give the universe zero spatial curvature, as deduced from the cosmic microwave 
background radiation Q]. 

Initial interpretation was to consider that in the universe there exists an important matter 



2 



component that, in its most simple description, has the characteristic of the cosmo logical 
constant A, i.e. a vacuum energy density which contributes to a large component of negative 
pressure, and thus accelerates rather than decelerates the expansion of the universe. An 
alternative interpretation is to consider quintessence dark energy, in the form of a scalar 
field with self interacting potentiality] . 

On the other hand, in recent years an important attention has received cosmological 
models in which bulk viscosity is considered. For example, it was shown that the introduction 
of this kind of viscosity into cosmo.ogica. models can avoid the big hang singularities^ fl 
13] and any contribution from particle production may be modelled as an effective bulk 
viscosity |8j]. The bulk viscosity arises typically in mixtures of different species (as in a 
radiative fluid) or of the same (but with different energies) fluids. The dissipation due to 
bulk viscosity converts kinetic energy of the particles into heat, and thus one expects it to 
reduce the effective pressure in an expanding fluid. This fact may play a crucial role in 
the inflationary era of the universe since it is interesting to know whether dissipative effects 
could be strong enough to do a large negative effective pressure leading to inflation. 

Many, probably most, of the inflationary cosmological models considered are of flat 
Friedmann-Robertson- Walker (FRW) type. In this kind of models the velocity gradients 
causing shear viscosity and temperature gradients leading to heat transport are absent. 
Thus, any dissipation in FRW universe may be modelled as a bulk viscosity within a ther- 
modynamical approach. These dissipative processes, that consider the bulk viscosityj^, lp| . 
are compatible with the homogeneity and isotropy assumptions for the universe, and may 
play an important role in the early universe, specially before nucleosynthesis . 

Curvature models have been studied with an important matter component whose equation 
of state is given by p = — p/3. Here, the universe expands at a constant speedj3]. On the 
other hand, a flat decelerating universe model has been described^!. Also, following a 
similar idea in accelerating models, the dark energy potential Q] and the location of the 
first Doppler peak in the Cosmic Microwave Background spectrum^] have been studied. 
In this way, we could consider models where the starting geometry were other than that 
corresponding to the critical geometry, but at low redshift these are indistinguishable from 
flat geometries. 

In this paper we want to study universe models in which we use the causal thermody- 
namical theory for processes out of equilibrium Q]. The stable and causal thermodynamics 
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of Israel and Stewart replaces satisfactorily the unstable and non-causal theory of Eckart 17] 
and Landau and Lifshitzfisl]. The main goal in this paper is to study not vanishing curva- 
ture causal inflationary universe models. We pretend to do this in Brans-Dicke (BD)-type 
theory Q]. Here the theory is characterized by a scalar field $ and a constant coupling 
function uj together with a scalar potential associated to From now on the subscript 
zero will represent the actual values. 

II. THE FIELD EQUATIONS 

We consider the effective action to be 

S = J d*x^rg{±-{$>R - ^{d^f + Vm + \{d,Qf + V(Q) + C M } (1) 

where $ is the BD scalar field, V A ($) is a scalar potential associated to the BD field, uq, is 
the BD parameter, Q the quintessence scalar field with associated potential V(Q), R the 
scalar curvature and Cm represents the matter contributions other than the Q component. 

We will obtain the field equations for a cosmological viscous model where the energy- 
momentum tensor of the fluid with a bulk viscosity is taken to be 

T^u = (p M + Vm + tyv^Vu + (p M + U)g lu/ , (2) 

with p M , p M the thermodynamical density, the pressure of the fluid, and II is the bulk viscous 
pressure, respectively. In the co-movil reference system we take v^v^ = — 1 and v^ 1 — 5q. 

The contribution to the energy-momentum tensor related to the quintessence scalar field, 
becomes given by: 

T^u = (p Q + Pq)u m u„ + g^Pq, (3) 
We can get the field equations of the universe from the action ( fTj) . to the FRW metric. 

2H + 3ff 2 + - 2 = ( I) - i($ + 2Hi>)+ 
a 2 2 \ $ / $ 
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(Pm - 3p M + p Q - 3P Q - 311) (6) 



2u; + 3 

where a(t) is the scale factor, and k is the curvature parameter that is equal to 1 (closed), 
(flat) or -1 (open). 

The matter conservation equation is given by 

p M + 3H{p M +p M + U) =0 (8) 

where the matter has a barotropic equation of state p M = (7 — l)p M where the parameter 7 
is in the ranger < 7 < 3. The energy density associated to the quintessence scalar field is 
given by, 

Pq = \Q 2 + V(Q) (9) 

and the pressure 

p q = \q 2 -v(Q). (10) 

This quintessence scalar field has a equation of state defined by p Q = w Q p Q , where w Q is 
the eq. of it parameter, and similar to the matter fluid, it satisfies the conservation equation 

p Q + 3H(p Q +p Q ) = 0. (11) 



III. CHARACTERISTICS OF THE MODEL 

In order to mimic a flat universe model, we impose the following conditions [l^. 

£*-£ + 5S v W = ° < 12) 



and 



$ 3 d (V(*)\ . 
pQ ~ 8tt(1 - 3^ Q ) d<t> V $ 2 )' 1 ' 
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In the following we will assume a power law potential for the BD field, V(<&) = Vo(J^)^, 
where Vq and (3 are constant. We will see that (3 becomes determined entirely by the equation 
state parameter uj q . With this assumption we obtain from the two latter eqs. 

k {2(3 -3(l + w Q )\ V /$ 



,2 



(14) 



6(1 -3w Q ) / $ l*o, 
where for V > and $ > 0, we obtain that (3 > |(1 + w Q ) if k > 0, (3 = |(1 + w Q ) if 
= and (3 < |(1+m; q ) if k < 0. Since we are interested in power law inflationary universe 
models, then we take the scale factor to be given by 

t^ N 



a{t) = a ( — J • (15) 



with N >1. Therefore, we obtain for the BD field 



2N 

$(t) =$o(r) (16) 



and Vo = 2^-3(14^) f^- With these solutions, the quintessence energy density results in 

2iV|3 

fl 6fc(^-2) $o/t\ — 

8vr(2/?-3(l+ WQ ))a 2 Uo7 1 J 

Since N > 1, then we need /3 < or (3 > 1 in order to satisfy the energy conservation 

equation. From the Q field eq.(JTTJ) this gives w Q = — !fz§> or equivalently (3 — 3 



Ifzf , or equivalently /3 = , 
which becomes determined in terms of the w Q parameter as was mentioned above. We 
should note that since — 1 < w Q < —1/3, then f3 < 0, therefore the possibility of satisfying 
f3 > 1 is not possible. This implies that the BD scalar potential is an inverse power law 
potential of the BD field. 

The field solutions together with the constraint equations yield to the following expression 
for the quintessence field: 



l+(JV-l)/3 



Q(t) = Qo[f] 1-/3 (is) 



where Q = J 8 -f 2 flitn +il jv - nffi w • From expressions © and together with the 



87r(2/3-3(l+u. Q ))(l+(A r -l)/3) 2 a 

equation of state for the quitessence field yields to 

2N/3 

V{Q) = V(Q ) (^-j (19) 

where V(Q ) = iQ^p-l^+wQ^ao 2 • Notice that, since /3 < 0, then we need just only to satisfy 
iV > 2 in order to get an usual quintessence scalar potential, i.e., an inverse power law scalar 
potential. 
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IV. TRANSPORT EQUATION FOR BULK VISCOSITY 

At this point we word like to describe dissipation due to bulk viscosity £, via Israel- 
Stewart theory. The bulk viscous pressure II is given by the transport equation (linear in 

n) 

rfl + n = -3H{ - irn Uh + T - - i - ^ V (20) 

where r is the relaxation time (which removes the problem of infinite propagation speeds) 
and T is the temperature of the fluid. In the non-causal formulation r = and then Eq. (J20|) 
has a simple form II = —3H!;. 

Following refs.0, [^(J we will take the thermodynamic quantities to be simple power 
functions of the matter density p M , 

i = ap M m ,T = pp M r y T = -^ = ap M m ~ 1 (21) 

Pm 

with a,p,m and r greater them zero. The expression for r is used as a simple procedure to 
ensure that the speed of viscous pulses does not exceed the speed of light. For an expanding 
cosmological model the constant m should be positive and we should satisfy 

r > H~\ (22) 



in order to have a proper physical behavior for £ and r jQL l21|. 

/,From the expansion law(|15|l and from eqs.(JHJ), (f2T)|) and (|2 1|) (with the standard relation 
for the temperature of a barotropic fluid in which r = (7 — 15/7) we obtain 

Pm + (3N+l)Pf + - P ]^ m p M + 3 -^ P 2 - m 
t a at 

-9iV 2 (l-l)^- (27 ~ 1) ^ = 0. (23) 
V 2) t 2 2 7 p M 1 ; 

This differential eq. may be solve for some specific values of the m parameter. Taking 
the following anzatsj2l| 

^ \ ( l-m ) 



"-=*U>) ■ (24) 

with m < 1 in order that p M decreases when t increase, we obtain from eq. (|23jl the constant 
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p° M becomes given by 



o = / g ) (25) 



(27)" 1 - 3iV(l - m) - 9iV 2 (l - 7 /2)(l - m) 2 \ 1/(1_m) 



(1 - 3^7(1 -m)) 

From this we can find a as a function of the parameters 7, N, m and p° . The pressure p M 
is obtained using the equation of state p M = (7 — l)p M . 

Now, from eq.(jHJ) and the barotropic eq. of state for the matter fluid we obtain for the 
bulk viscous pressure 

n = -|^-7P M , (26) 



and substituting eq.(j24j) in this latter eq. we get 

We see that IT < if 7 > l/(3iV(l - m)), n > if 7 < l/(3iV(l - m)) and n = if 
m = 1 - 1/(3A^7). 

In order to satisfy the field eqs. of motion the parameters should follow the equality 
= 2 — with (3 expressed in terms of w Q , thus m becomes a function of the parameter 
N and w Q , 

11-3 

m = 1 — 28 

21-P-N y 1 

where 3 is a function of the w Q parameter. 

At local equilibrium, the entropy satisfy the Gibbs equation 

TdS = (p M + p M )d (-)+ -dp M , (29) 
\n J n 

where n is the number density and satisfy the conservation equation 

h + 3Hn = 0, (30) 

from which we get 

t ,-3N 



n(t) = no ( - J (31) 



or equivalently n(t)a 3 (t) = const. 
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Using (J29|) and (J31|) we obtain the well known evolution equation for the entropy (neglect- 
ing the heat flux and the shear viscosity) 

thus from eqs. flHJ, dZH>, (Ei, (E3) andfjSH) we get 

For proper physical behaviour we must satisfy S > 0, which implies the condition m < 
1 — l/j(3N — 1). This inequality together with eq.([28|) give an inequality for 7 as function 
of the parameters N and w Q , 

7>2 _ (i-p-N) 



;i - ;3)(3/V - 1) 
Thus, this parameter becomes bounded from below. 

The total entropy E in a comoving volume is defined by means E = Sna 3 (t). Then, by 
taking Eq. (j31j) and()32|) we may write the growth of total nondimensional comoving entropy 
over a proper time interval tf and ti asjjj 

Er^-^'^ (34) 

Now taking ti and tf as the beginning and the exit time for the inflation period of the 
universe respectively and from eqs.(JT3J), (J3TJ), (j2U), (}2"Tj) . and (|3*3|) we obtain for the increase 
the total nondimensional entropy in the comoving volume a 3 (t) the following expression 

L,f-L,i — 



k B \i 



(35) 



We take the values for the beginning and ending time of inflation as ti ~ 10~ 35 s and 
tf ~ 10 -32 s, respectively. In the following numerical calculation we take the reference 
time to equal to tf. Considering that the universe at the end of the period inflation exits 
to the radiation era, we can constraint some of the constants of integration of the above 
formulae. Effectively, it is known that the temperature of the universe at the beginning of 
the radiation era is approximately T ~ 10 14 GeV ~ 1.16 x 10 27 K. Thus, the temperature 
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at the end of inflation must be Tf = 1.16 x 10 27 K. On the other hand, we know that in 
this period p = a r T 4 , where a r = 15c3 ^ ~ 7.56 x 10~ 15 J m^K~ A . From here we conclude 
that at the end of inflation any period (or at the beginning of the radiation era) we have 
p f « 10 93 Jm~ 3 . This implies that r = 1/4 for the exponent in the temperature T, i.e. 
7 = 4/3, as we can see from eq. [2T] 

Thus the typical values for the inflationary period are[9j: 

U « 10" 35 s; t f w 10~ 32 s; a { « cfc<, 
T/ « 10 27 fT, 7 = 4/3, p « x 10 93 J/m 3 . (36) 

The e-folding parameter Z = ln[a(tf)/a(ti)\ for the power law inflation (J15j) takes the form 



Z = nln\^-J. (37) 

It is well known that for solving most the problems of the standard model in cosmology we 
must have Z « 60 - 70. Thus from lj2I) . l|55|) . (f36l andflHH) we have 

^ e3Z ^ 3/ ^l-(10-3)--^T), (38) 



3fc B /i 

where we have used the relation a 3 = a 3 e 3Z , following from the definition of Z, eqJ37l We 
see from Eq. (}3*Hj) that if the inequality 3A^7(1 — m) » 1 is satisfied (in this case Ej >> E$) 
we obtain the accepted value for the total entropy in the observable universe jl^] 

E^xlO 88 . (39) 

Then this model can account for the generally accepted entropy production via causal dis- 
sipative inflation, without reheating. 

Finally, for complete, we will see what happen when m — 1. For this value and considering 
eq.(J2HJ), which gives the solution j^jj 

, . v -7(3iV+2) 

Pu = P W t,a (l) , (40) 

where N 2 = 27/9 and p° M is an arbitrary constant. Unfortunately this solution implies a 
decreasing entropy 

which certainly we should neglect in order to finish with a reasonable physical model. 
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V. CONCLUSIONS 

To conclude, on the light of the above results we may say that there exist an intimate 
relation between causal dissipative inflation and reheating. We expect to come back to this 
point of research in the near future, where extensions to more general universe models will 
be applied. 
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